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We study the correlation between the thickness of the neutron skin in finite nuclei and the nuclear 
symmetry energy for isotopic chains of even-even Ni, Sn, and Pb nuclei in the framework of the 
deformed self-consistent mean-field Skyrme HF-I-BCS method. The symmetry energy, the neutron 
pressure and the asymmetric compressibility in finite nuclei are calculated within the coherent 
density fluctuation model using the symmetry energy as a function of density within the Brueckner 
energy-density functional. The mass dependence of the nuclear symmetry energy and the neutron 
skin thickness are also studied together with the role of the neutron-proton asymmetry. A correlation 
between the parameters of the equation of state (symmetry energy and its density slope) and the 
neutron skin is suggested in the isotopic chains of Ni, Sn, and Pb nuclei. 
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I. INTRODUCTION 

The nuclear symmetry energy is a quantity of crucial 
importance in different areas of nuclear physics, includ- 
ing structure of ground-state nuclei [l|-|3|, dynamics of 
heavy- ion reactions j^-Q : physics of giant collective exci- 
tations Q and physics of neutron stars [sl-fToj. Recently, 
the interest in the symmetry energy has been stirred up 
by novel astrophysical observations and by the availabil- 
ity of exotic beams in accelerators that provide additional 
information to the standard nuclear asymmetry studies 
based on stable nuclei. Particularly important in the 
different areas, and similarly uncertain, is the density 
dependence of the symmetry energy in uniform matter. 
As can be seen e.g., in Refs. |lll4l5| . an increasing wide 
range of theoretical conclusions are being proposed on 
that density dependence, as well as on some associated 
nuclear characteristics. In the last years, the tempera- 
ture dependence of single-particle properties in nuclear 
and neutron matter was also broadly investigated (e.g., 

Refs. 

Measurements of nuclear masses, densities, and col- 
lective excitations have allowed to resolve some of the 
basic features of the equation of state (EOS) of nuclear 
matter. However, the asymmetry properties of the EOS 
due to different neutron and proton numbers remain 
more elusive to date, and the study of the isospin de- 
pendent properties of asymmetric nuclear matter |18l422l | 
and the density dependence of the nuclear symmetry en- 
ergy remains a prime objective. The new radioactive 
ion beam facilities at GSR (Ghina), FAIR (Germany), 
RIKEN (Japan), SPIRAL2/GANIL (France) and the up- 
coming FRIB (USA) will provide the possibility of ex- 
ploring the properties of nuclear matter and nuclei under 
the extreme condition of large isospin asymmetry. 

Experimentally, the symmetry energy is not a directly 
measurable quantity and has to be extracted indirectly 



from observables that are related to it (see, for exam- 
ple, the recent review [11]). The neutron skin thickness 
of nuclei is a sensitive probe of the nuclear symmetry 
energy, although its precise measurement is difficult to 
obtain. At present neutron skin thicknesses are derived 
from pygmy dipole resonances measurements [23|, data 
from antiprotonic atoms [25i] and other methods for its 
extraction like reactions and giant resonances. This al- 
lows one to constrain the parameters describing the nu- 
clear symmetry energy. 

Considering nuclei from the isotopic chains of Sn and 
Pb, Warda et al. [1^ studied theoretically the bulk and 
the surface nature of the formation of the neutron skin, a 
concept that can be applied when analyzing the experi- 
mental data. The same authors indicated in Ref. j27j] the 
role of the stiffness of the nuclear symmetry energy on 
the origin of the neutron skin thickness of ^°®Pb, the lat- 
ter being decomposed into bulk and surface components. 
Also, Danielewicz [2^ has demonstrated that the ratio 
of the volume symmetry energy to the surface symmetry 
energy is closely related to the neutron skin thickness. 
The correlation between bulk and surface symmetry en- 
ergy was further discussed in Refs. [1, [29[ - l32j . Recent 
calculations of thermal nuclear properties showed that 
the surface symmetry energy term is more sensitive to 
temperature than the volume energy term [s^ . 

The neutron skin thickness, generally defined as the 
difference between neutron and proton rms radii in the 
atomic nucleus, is closely correlated with the density de- 
pendence of the nuclear symmetry energy and with the 
equation of state of pure neutron matter. Moreover, it 
has been shown that the neutron skin thickness in heavy 
nuclei, like ^°*Pb, calculated in mean-field models with 
either nonrelativistic or relativistic effective nuclear in- 
teractions, displays a linear correlation with the slope of 
the neutron EOS obtained with the same interactions at 
a neutron density p 0.10 fm~^ [sl, [s^. The statistical 
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analysis performed in Ref. 13 6| clearly shows the strong 
correlation that exists between the neutron form factor of 
^°^Pb and isovector indicators such as neutron skins and 
radii in neutron-rich nuclei, dipole polarizability, and nu- 
clear matter properties (symmetry energy and pressure). 

The symmetry energy of finite nuclei at saturation den- 
sity is often extracted by fitting ground state masses with 
various versions of the liquid-drop mass (LDM) formula 
within liquid-drop models [stI - IsqI I . It has been also stud- 
ied in the random phase approximation based on the 
Hartree-Fock (HF) approach j4^] or effective relativistic 
Lagrangians with density-dependent meson-nucleon ver- 
tex functions (4]| . energy density functionals of Skyrme 
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force 42j-j44| as well as relativistic nucleon-nucleon inter- 



46| . In the present work the symmetry energy 
will be studied in a wide range of finite nuclei on the ba- 
sis of the Brueckner energy-density functional for nuclear 
matter [47, 4^ and using the coherent density fiuctuation 
model (CDFM) (e.g., Refs. [ii,!!^). The latter is a natu- 
ral extension of the Fermi gas model based on the gener- 
ator coordinate method (sol . [5l| and includes long-range 
correlations of collective type. The CDFM has been ap- 
plied to calculate the energies, the density distributions 
and rms radii of the ground and first monopole states in 
■*IIe, and '^^Ca nuclei 52]. An analysis of breath- 
ing monopole states whose description is related mainly 
to general characteristics of atomic nuclei and weakly to 
the peculiarities of nuclear structure in the framework 
of the CDFM has been performed and also results for 
the incompressib ility of finite nuclei have been reported 
[Hil. In Refs. Ilii] the CDFM has been employed to 
calculate the scaling function in nuclei using the relativis- 
tic Fermi gas scaling function and the result has been ap- 
plied to lepton scattering processes [53 - [60| . The CDFM 
analyses became useful to obtain information about the 
role of the nucleon momentum and density distributions 
for the explanation of superscaling in lepton-nucleus scat- 
tering [ssl, H^. The CDFM scaling function has been 
used to predict cross sections for several processes: inclu- 
sive electron scattering in the quasielastic and A regions 
[stI [ssf and neutrino (antineutrino) scattering both for 
charge-changing [H, [g^ and for neutral-current [H, |6^ 
processes. Recently, the CDFM was applied to study 
the scaling function and its connection with the spectral 
fimction and the momentum distribution 61]. 

In our previous study [6^ . we analyzed within a de- 
formed HF-I-BCS approach with Skyrme-type density- 
dependent effective interactions the existence of a skin in 
nuclei far from the stability line by testing various defi- 
nitions previously proposed (63l |64| . We have found that 
the neutron skin determined by the difference between 
neutron and proton radii using diffraction parameters 
defined in the Helm model [63] exhibits a more smooth, 
gradual increase with the neutron excess than the one ob- 
tained by means of definition from Ref. 6J]. As a result, 
the absolute size of the skin ranges from 0.4 fm when the 
Helm model is used up to almost 1 fm for the heaviest 
isotopes using the criteria from Ref. (64| . The same Helm 



model definitions of the neutron skin thickness were stud- 
ied in Ref. ^] in spherical self-consistent Hartree-Fock- 
Bogoliubov (HFB) calculations with finite-range Gogny 
forces. We would like to emphasize the consistency of 
the results of both approaches for the case of Ni and Sn 
isotopes [g^, HE]. Some particular investigations of the 
neutron skin thickness of neutron-rich nuclei and its con- 
nections with the symmetry energy have been performed 
in the droplet model with surface width dependence [g^ 
and in studies based on the Landau-Migdal approxima- 
tion [13]. 

In the present work we investigate the relation between 
the neutron skin thickness and some nuclear matter prop- 
erties in finite nuclei, such as the symmetry energy at 
the saturation point, symmetry pressure (proportional to 
the slope of the bulk symmetry energy), and asymmet- 
ric compressibility, considering nuclei in given isotopic 
chains and within a certain theoretical approach. In ad- 
dition to various linear relations between several quanti- 
ties in bulk matter and for a given nucleus that have been 
observed and tested within different theoretical methods 
(e.g. nonrelativistic calculations with different Skyrme 
parameter sets and relativistic models), we are looking 
forward to establish a possible correlation between the 
skin thickness and these quantities and to clarify to what 
extent this correlation is appropriate for a given isotopic 
chain. As in our previous papers [6^ [68{ , we choose some 
medium and heavy Ni and Sn isotopes, because, first, 
they are primary candidates to be measured at the up- 
coming experimental facilities, and second, for them sev- 
eral predictions have been made concerning the nuclear 
skin emergence. In addition, we present some results for 
a chain of Pb isotopes being inspired by the significant in- 
terest (in both experiment [6^ and theory 7.0] ) to study, 
in particular, the neutron distribution of ^"*Pb and its 
rms radius. The densities of these nuclei were calculated 
within a deformed HF-I-BCS approach with Skyrme-type 
density-dependent effective interactions [7l|, [72| . The re- 
sults obtained in Ref. [g^] demonstrated the ability of 
our microscopic theoretical method to predict the nu- 
clear skin in exotic nuclei. As already mentioned, for a 
link between nuclear matter and finite nuclei we adopt 
the energy-density functional of Brueckner et al. [13, 13 
that has been able to reproduce the binding energies and 
rms radii of light, medium, and heavy nuclei. The intrin- 
sic properties of these nuclei are finally calculated in the 
CDFM framework that provides an extension to realis- 
tic finite nuclear systems. The evolution of the symmetry 
energy as we increase the number of neutrons is also stud- 
ied. Additionally to Ref. [62|, we pay more attention to 
the role of the neutron-proton asymmetry on the neutron 
skin formation. 

The structure of this article is the following. In Sec. II 
we present the common definitions of symmetry energy 
and properties of nuclear matter which characterize its 
density dependence around normal nuclear matter den- 
sity. Section HI contains the CDFM formalism that pro- 
vides a way to calculate the intrinsic quantities in finite 



3 



nuclei. There we give also the basic expressions for the 
model density distributions and nuclear radii obtained in 
the deformed HF+BCS method. The numerical results 
and discussions are presented in Sec. IV. We draw the 
main conclusions of this study in Sec. V. 



II. THE KEY EOS PARAMETERS IN 
NUCLEAR MATTER 

The symmetry energy S{p) is defined by a Taylor series 
expansion of the energy per particle for nuclear matter 
(NM) in terms of the isospin asymmetry 5 = (/o„ — pp)/ p 



E{p, S) = E{p, 0) + S{p)S^ + 0(3^) + 



(1) 



where p = pn + Pp is the baryon density with p„ and pp 
denoting the neutron and proton densities, respectively 
(see, e.g. [13, [z3)- Odd powers of 5 are forbidden by the 
isospin symmetry and the terms proportional to (5^ and 
higher orders are found to be negligible. 

Near the saturation density po the energy of isospin- 
symmetric matter, E{p,0), and the symmetry energy, 
S{p), can be expanded as 



E{p,0)=Eo + -^ip- Po)^ + 



and 



1 d^E{p,6) 



5=0 



(2) 



a4 + — (p - pq) + T^(P - po) + 
Po l^Po 



(3) 



The parameter is the symmetry energy at equilibrium 
(p = Po). The pressure -p^^^ 



NM _ 2 9S 



(4) 



a factor of three) which is important for structure calcu- 
lations. In practice, predictions for the symmetry energy 
vary substantially: e.g., 04 = S'(po) = 28 — 38 MeV. An 
empirical value of 04 k, 29 MeV has been extracted with 
reasonable accuracy from finite nuclei by fitting ground- 
state energies using the generalized Weizsacker mass for- 
mula ii: 



E{N, Z) 



Ey + Es + Ea + Ec + E,nic 



+ ac 



-ay A + asA"^/^ 



A 



E„ 



(7) 



where N and Z are the neutron and proton numbers, re- 
spectively, and A = N+Z is the mass number. In Eq. ([7]) 
ay, as, fla and ac are the volume, surface, symmetry and 
Coulomb coefficients. The microscopic energy Emic con- 
tains shell- and pairing-energy corrections. It is worth 
to mention that the volume component of the symmetry 
energy coefficient Oa (or volume-symmetry coefficient in 
the empirical mass formula ([7])) has to be attributed to 
the parameter 04 in Eq. ([3]). 



III. SYMMETRY ENERGY PARAMETERS OF 
FINITE NUCLEI IN CDFM 

The CDFM was suggested and developed in Refs. [i^, 
[soj . The model is related to the delta- function limit of 
the generator coordinate method (sol . [sij . It is shown in 
the model that the one-body density matrix of the nu- 
cleus p(r,r') can be written as a coherent superposition 
of the one-body density matrices (OBDM) for spheri- 
cal "pieces" of nuclear matter (so-called "fluctons") with 
densities 



Pxir) = po(a;)e(a; - |r|), 



(8) 



and the curvature AK 



NM 



AK 



NM 



9Po 



dp^ 



(5) 



P=Pq 



of the nuclear symmetry energy at po govern its density 
dependence and thus provide important information on 
the properties of the nuclear symmetry energy at both 
high and low densities. The widely used "slope" param- 
eter L^^^ is related to the pressure p^*^ [Eq. dU] by 



j^NM ^ "^Po 



NM 



Po 



(6) 



We remark that our present knowledge of these ba- 
sic properties of the symmetry term around saturation 
is still very poor (see the analysis in Ref. 7J| and refer- 
ences therein). In particular we note the uncertainty of 
the symmetry pressure at po (sometimes it can vary by 



where 



Po{x) 



iA 



(9) 



The generator coordinate x is the radius of a sphere con- 
taining Fermi gas of all A nucleons uniformly distributed 
in it. It is appropriate to use for such a system OBDM 
of the form; 



Px{r,r') = 3po(x) 

X e( X- 



ji(fcF(x)|r - r'l) 
(fcF(x)lr-r'l) 
|r + r' 



(10) 

where ji is the first-order spherical Bessel function and 

'3^2 X 1/3 



kpix) = ( — Po(a;) 



a 

X 



(11) 
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with 



where 



V 8 



(12) 



is the Fermi momentum of such a formation. Then the 
OBDM for the finite nuclear system can be written as a 
superposition of the OBDM's from Eq. (ITU|) : 

/>oo 

p{r,r')= / dx|^(x)|Vx(r,r'). (13) 
Jo 

In the CDFM the Wigner distribution function which 
corresponds to the OBDM from Eq. (|13l) is: 



/>oo 

W{r,k)= dx\J^{x)fW,{r,k), (14) 
Jo 



where 



VK,(r,k) 



(2^) 



.e(x-|r|)e(fc;^(a:)-|k|). (15) 



Correspondingly to VF(r,k) from Eq. (|14l) . the density 
p(r) in the CDFM is expressed by means of the same 
weight function |J^(a:)p: 



p(r) 



dkW{r,k) 

DO 



3A 

Attx^ 



normalized to the mass number: 



J p{r)dr = A. 



Q{x-\r\) (16) 



(17) 



If one takes the delta-function approximation to the Hill- 
Wheeler integral equation in the generator coordinate 
method one gets a differential equation for the weight 
function J-{x) [13, [5l|. Instead of solving this differ- 
ential equation we adopt a convenient approach to the 
weight function | J"(a:)p proposed in Refs. [4^,[5^. In the 
case of monotonically decreasing local densities {i.e. for 
dp{r)/dr < 0), the latter can be obtained by means of a 
known density distribution p{r) for a given nucleus [from 
Eq. dm)]: 



\:f(x) 



1 



dp{r) 



po{x) dr 

The normalization of the weight function is: 



dx\J-{x) 



1. 



(18) 



(19) 



Considering the pieces of nuclear matter with density 
Po{x) (HJ one can use for the matrix element V{x) of the 
nuclear Hamiltonian the corresponding nuclear matter 
energy from the method of Brueckner et al. [U 13 ■ 
this energy-density method the expression for V{x) reads 



Vix)=AVo{x) + Vc-Vco, 



(20) 



Voix) = 37.53[(l + J)'^/3^(l-(5)5/3]^2/3^^^ 

+ bipo{x) + h2py^{x) + h'ipy^{x) 

+ 5^[h^po{x) + hpl'\x) + KpI'\x)] (21) 

with 

bi = -741.28, 62 = 1179.89, h = -467.54, 

64 = 148.26, &5 = 372.84, fog = -769.57. (22) 

Vo(x) in Eq. (j20l) corresponds to the energy per nu- 
cleon in nuclear matter (in MeV) with the account for 
the neutron-proton asymmetry. Vc is the Coulomb en- 
ergy of protons in a "flucton": 



Vc 



5 x 



and the Coulomb exchange energy is: 

Vco 0.7386Ze2(3Z/47rx3)i/^ 



(23) 



(24) 



Thus, in the Brueckner EOS [Eq. the poten- 

tial symmetry energy turns out to be proportional to 5^ . 
Only in the kinetic energy the dependence on 5 is more 
complicated. Substituting Vo{x) in Eq. ^ and taking 
the second derivative, the symmetry energy S^^^{x) of 
the nuclear matter with density po{x) (the coefficient 04 
in Ec^. ^) can be obtained: 



:,NM I 



{x) = 41.7p^/'(.T)+&4Po(x) 

+ hpl'\x) + Kpl'\x). (25) 

The corresponding analytical expressions for the pressure 
Pq^^{x) and asymmetric compressibility AK^^''^{x) of 
such a system in the Brueckner theory have the form: 



NM 

Po 



(x) 



= 27 
4 

+ 



^o^^{x) + biplix) 
^hpl^^ix) + -b&pl'^^x) 



(26) 



and 



Aif^*^(x) = -&Z.Apl"'{x)+Ab^py\x) 



mQpl'\x). 



(27) 



Our basic assumption within the CDFM is that the 
symmetry energy, the slope and the curvature for finite 
nuclei can be defined weighting these quantities for nu- 
clear matter (with a given density po{x) ([9])) by means 
of the weight function |J"(2;)|^. Thus, in the CDFM they 
will be an infinite superposition of the corresponding nu- 
clear matter quantities. Following the CDFM scheme, 
the symmetry energy, the slope and the curvature have 
the following forms: 



dx\F{x)YS'^^\x) 



(28) 
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po = I dx\J^ix)\^p^''{x), 
/o 



AK^ dx\J^{x)\^AK^^\x). 
Jo 



(29) 



(30) 



We note that in the limit case when p{r) = poQ{R—r) (po 
being the saturation density of symmetric nuclear mat- 
ter) and |J^(x)p becomes a delta- function (see Eq. ^T8\i ). 
Eq. (Eg) reduces to 5^*^(/9o) = 04- 

The results that will be shown and discussed 
in the next section have been obtained from self- 
consistent deformed Hartree-Fock calculations with 
density-dependent Skyrme interactions [7l| and account- 
ing for pairing correlations. Pairing between like nucleons 
has been included by solving the BCS equations at each 
iteration with a fixed pairing strength that reproduces 
the odd-even experimental mass differences. 

We consider in this paper the Skyrme SLy4 [tI], Sk3 
[tg} and SG2 [t^] parametrizations, namely because they 
are among the most extensively used Skyrme forces and 
have been already used in our previous paper [6^ . In 
addition, we probe a recent LNS parameter set from 
Ref. [t^, derived in the framework of the Bruekner- 
Hartree-Fock approximation that has been shown to work 
reasonably well for describing finite nuclei in the HE ap- 
proximation. 

The spin-independent proton and neutron densities are 



i ne spm- mdep ' 
given bv [gIIt^ 

p{R) =^ p{r, z)^^ '2vfpi{r, z) 



(31) 



in terms of the occupation probabilities v1 resulting from 
the BCS equations and the single-particle densities pi 



pdR) = P^{r. z) = m{r, z)\' + |<i>r(r, z)\' , (32) 



with 



$±(r,z) 



27r 



X ^s,±i/2 -^A.AT V'l ir) V-™, (z) (33) 

a 

and a = {ur, Uz, A, S}. The functions ip^^ (r) and tpn^ (z) 
entering Eq. (j33p are defined in terms of Laguerre and 
Hermite polynomials 



(rir + A)! 



(3^V2r^^f'e-^/'LUv), (34) 



(-) = ^/t^^ ^'J' Hn. (C) , (35) 



with 



e = z/3. , i^^r^Pl. (36) 



The normalization of the densities is given by 

p{R)dR = X, (37) 



with X = Z, N for protons and neutrons, respectively. 

The mean square radii for protons and neutrons are 
defined as 



^ _ S R^p,,AR)dR 
JppAR)dR 

and the rms radii are given by 

'p,n — 'p^n ^ 



(38) 



(39) 



Then the neutron skin thickness is usually characterized 
by the difference of neutron and proton rms radii: 



AR =< rl >i/2 _ < ^2 ^1/2 



IV. RESULTS AND DISCUSSION 



(40) 



As the main emphasis of the present study is to inspect 
the correlation of the neutron skin thickness AR of nuclei 
in a given isotopic chain with the s, po and AK param- 
eters extracted from the density dependence of the sym- 
metry energy around saturation density, we show first in 
Fig. [H the results for Ni isotopes. The symmetry energy, 
the pressure and the asymmetric compressibility are cal- 
culated within the CDEM according to Eqs. (pg|) - ((5g)) 
by using the weight functions ((T8| calculated from the 
self-consistent densities in Eq. (|5T|) . The differences be- 
tween the neutron and proton rms radii of these isotopes 
[Eq. pO)) ] are obtained from HF+BCS calculations using 
four different Skyrme forces, SLy4, SG2, Sk3 and LNS. 
It is seen from Fig. [TJa) that there exists an approximate 
linear correlation between AR and s for the even-even 
Ni isotopes with A = 74 — 84. We observe a smooth 
growth of the symmetry energy till the double-magic nu- 
cleus ^^Ni {N = 50) and then a linear decrease of s while 
the neutron skin thickness of the isotopes increases. This 
behavior is valid for all Skyrme parametrizations used in 
the calculations, in particular, the average slope of AR 
for various forces is almost the same. The LNS force 
yields larger values of s comparing to the other three 
Skyrme interactions. In this case the small deviation can 
be attributed to the fact that the LNS force has not been 
fitted to finite nuclei and therefore, one cannot expect a 
good quantitative description at the same level as purely 
phenomenological Skyrme forces [78j. As a consequence, 
the neutron skin thickness calculated with LNS force has 
larger size with respect to the other three forces whose 
results for AR are comparable with each other. Never- 
theless, it is worth to compare its predictions not only for 
ground-state nuclear properties but also for EOS charac- 
teristics of finite systems with those of other commonly 
used Skyrme forces. 
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FIG. 1: (Color online) HF+BCS neutron skin thicknesses AR for Ni isotopes as a function of the symmetry energy s (a), 
pressure po (b), and asymmetric compressibility AK (c) calculated with SLy4, SG2, Sk3, and LNS forces. 



We also find a similar approximate linear correlation 
for Ni isotopes between AR and po [Fig. [Hb)] and less 
strong correlation between AR and AK [Fig. [IJc)]. As 
in the symmetry energy case, the behavior of the curves 
drawn in these plots shows the same tendency, namely 
the inflexion point transition at the double-magic ^^Ni 
nucleus. We would like to note that the predictions for 
the difference AR between the rms radii of neutrons and 
protons with Skyrme forces obtained in Ref. [63| exhib- 
ited a steep change at the same place in which the num- 
ber of neutrons starts to increase in the chain of nickel 
isotopes. In addition, the calculated results for the two- 
neutron separation energies and neutron and matter rms 
radii of the even Ni isotopes obtained in the relativistic 
Hartree-Bogoliubov framework [t^ showed a quite strong 
kink at iV = 50. Also one can see from Figs. Hfb) and 
[TJc) that the calculated values for po and AK are lower 
in the case of LNS force than for the other three Skyrme 
parameter sets. 

The results presented in Fig. [TJa) are illustrated in 
more details in Fig. [51 where the evolution of both neu- 
tron skin thickness and symmetry energy for Ni isotopes 
calculated with SLy4 force is given simultaneously on the 
same figure. For this chain our predictions for AR are 
extended to more neutron-rich isotopes {A = 84) com- 
paring to the range of isotopes considered in Ref. [6^ . 
The results for the symmetry energy obtained from the 
four different Skyrme parametrizations are compared in 
Fig. 131 It is seen that the trend of s when the mass num- 
ber increases is preserved for all Skyrme forces but its 
values vary being very similar for SG2 and Sk3 forces. 
However, the magnitude of the symmetry energy values 
slowly changes and turns out to be approximately in the 
range of 27-29 MeV. The nuclear matter calculations re- 
lying on realistic nucleon-nucleon interactions appear to 
yield volume symmetry energy values within the same 
range of 27-29 MeV depending on the interaction ^8d\. 

As is known, radioactive ion beams enable us to study 
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FIG. 2: HF-I-BCS neutron skin thicknesses AR (solid line) 
and symmetry energies s (dashed line) for Ni isotopes calcu- 
lated with SLy4 force. 




FIG. 3: (Color online) The symmetry energies s for Ni iso- 
topes calculated with SLy4, SG2, SkS, and LNS forces. 
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radii of unstable nuclei over a wide range of the relative 
neutron excess (or relative neutron-proton asymmetry) 
/ ~ {N — Z) I A. The role of the global asymmetry on the 
neutron skin thickness in the Ni isotopic chain is shown in 
Fig- m for the case of SLy4 force by presenting results for 
the ratio Ai?/7 as a function of the mass number A. It 
can be seen from Fig.|4]that due to the gradual increase of 
the skin thickness almost no dependence on / is observed 
for isotopes with A < 78, while for the heavier ones the 
ratio Ai?// starts to increase rapidly with A. The latter 
takes place namely in the region where a pronounced 
neutron skin in Ni isotopes can be expected. We also 
give in Fig. |4] estimations for this ratio using the extended 
LDM formula for the binding energy of the nucleus [H, 
l30l | (with neglected Coulomb contribution): 



LDM 



A^R 



6NZ{1 + A^/^/ysY 



(41) 



In Eq. (|4T|l ys = Sy/Ss is the ratio between the vol- 
ume and surface synimetry energies and its value is taken 
to be ys = 1.7 js^]. For R we use R ^ rpA'^/^ with 
ro = 1.2 fm. As pointed out by Danielewicz [28], Eq. (|4T|) 
is valid for the difference of sharp-sphere radii, but it is 
generalized for the case of rms radii that must be used. 
This creates an additional Coulomb correction j2^ pre- 
serving unchanged the relation between AR and /. We 
see that the difference between the neutron and proton 
radii is linear in the global asymmetry in the absence of 
the Coulomb contribution (which, in principle, cannot be 
neglected) and measures the symmetry coefficient ratio. 
Thus, Fig. |4] displays how much the mean-field results 
differ from the LDM ones following the general obser- 
vation that the neutron skin depends roughly linearly on 
the asymmetry of the nucleus. Although the Coulomb ef- 
fects result in a reduction of the neutron skin for iV > Z 
nuclei (28l . [30| , it is important to note the direct relation 
between the skin thickness AR and the symmetry energy 
parameters {Sv,Ss) given by Eq. (|4T1) . 

In this aspect, here we would like to mention that the 
finite nucleus symmetry energy within the extended LDM 
(which accounts for the surface effects) is given by [28i - [30| 



aaiA) 



l + t/,A-i/3- 



(42) 



If one takes the same realistic values of Sy = 27 MeV and 
Vs = 1-7 [IqI which we use when dealing with Eq. (|4T|). 
it turns out that the values of are in the interval 19- 
21 MeV for the nuclei ranging between ^*^Ni and 2°**Pb. 
These values are significantly smaller than the symmetry 
energy values s derived in CDFM from Eq. ((28)) . As far 
as one of the main goals of our work is to calculate the 
nuclear symmetry energy of finite nuclei and to look for 
its relationship with the neutron skin thickness in these 
nuclei from a given isotopic chain, we would like to em- 
phasize that the results obtained by using Eq. (1^51) repre- 
sent estimations of the symmetry energy for finite nuclei 
obtained on the basis of nuclear matter EOS. A way to 




FIG. 4: The ratio AR/I for Ni isotopes calculated with SLy4 
force (solid line) and obtained from Eq. (|4ip with LDM [2a . 
[30I ] (dashed line). 



estimate the symmetry energy in finite nuclei within the 
mean-field framework has been recently proposed in [8l| , 
where Samaddar et al. studied the density and excita- 
tion energy dependence of symmetry energy and symme- 
try free energy of hot nuclei by using the local density 
approximation. 

In addition, it is important to note that by means of 
Eq. (j28p we obtain the symmetry energy in finite nu- 
clei including both bulk and surface contributions on the 
base of the Brueckner EOS. This is related to the us- 
age in the method of the weight function |J^(a;)p which 
is computed by taking the derivative of the density dis- 
tribution of a given finite nucleus [Eq. p^ ]. Thus, in- 
formation from both parts (central and surface) of the 
nuclear density is included in |J^(a;)|^ and, consequently, 
in the total symmetry energy. In this way, starting from 
the nuclear matter symmetry energy, we calculate by our 
method the total symmetry energy for the given finite 
nucleus. The use of different types of Skyrme effective 
forces leads to density distributions of protons and neu- 
trons which slightly modify the weight function |J^(a;)|^. 
Increasing the number of neutrons at fixed proton num- 
ber, Eq. (|28p provides values for the symmetry energies s 
of considered isotopes that are found to change smoothly 
within a given chain. 

As has been recently discussed (see, for instance, 
Refs. [2^, HO, 113); the relative contributions of the vol- 
ume symmetry and surface symmetry energies can be 
disentangled using the extended LDM formula for the 
binding energy of the nucleus. A comparison between 
the symmetry energy values obtained in [s^l (and listed 
in Table I of this reference) for SLy4 parametrization 
and those calculated in the present paper with the same 
force can be made. For example, the values of the to- 
tal symmetry energy for the case of SLy4 extracted from 
leptodermous expansion in Ref. [s^ and Eq. (j42p vary 
between 19 and 23 MeV. Our values of the total sym- 
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metry energy (deduced with from calculations of 

the density corresponding to SLy4 force [Eq. ([T^]) lie in 
the range of 26.8-27.8 MeV for Ni, 27.3-28.8 MeV for 
Sn, and 28.7-29.1 MeV for Pb isotopes. In our opinion, 
the differences of the compared results are due to the 
peculiarities of the two methods. As noted in the previ- 
ous paragraph, in Eq. (|28p we use obtained from 
both parts (central and surface) of the nuclear density, 
weighting by it the total symmetry energy for nuclear 
matter S^^{x). Second, in the present work S^^^{x) is 
taken from the Brueckner theory, while in (s^ (using the 
results of [sS]) the bulk symmetry energy in the consid- 
ered case of SLy4 parametrization is used. In principle, 
other functionals apart from the Brueckner one reported 
here can be applied within our method to calculate the 
total symmetry energy for finite nuclei, but as described 
above, the method gives it unseparated into bulk and sur- 
face contributions. 

In Fig. [5] we display the theoretical neutron skin thick- 
ness AR of nuclei from a Pb isotopic chain against the 
parameters of interest s, po and AK. Similarly to the 
three panels for Ni isotopes presented in Fig. [T] the pre- 
dicted correlations manifest the same linear dependence. 
However, in this case the kink observed at double magic 
208p|-, nucleus (/ ~ 0.21) is much less pronounced that 
it can be seen at ^^Ni isotope (/ « 0.28) and it does 
not change its direction. The value of AR for ^"^Pb 
deduced from present HF-I-BCS calculations with SLy4 
force is lower than the predicted thickness by the Skyrme 
Hartree-Fock mod el i 42il and by the extended relativistic 
mean- field model [4^, but fits well the values calculated 
with self-consistent densities of several nuclear mean-field 
models (see Table I in Ref. 01). The po and AK val- 
ues for ^"^Pb are in a good agreement with those from 
Ref. El. 

The analysis of the correlation between the neutron 
skin thickness and some macroscopic nuclear matter 
properties in finite nuclei is continued by showing the 
results for a chain of Sn isotopes. This is done in 
Fig. H where the resuhs obtained with SLy4, SG2, Sk3, 
and LNS Skyrme forces are presented for isotopes with 
A = 124 — 152. Similarly to the case of Ni isotopes with 
transition at specific shell closure, we observe a smooth 
growth of the symmetry energy till the double-magic nu- 
cleus ^^^Sn (iV = 82) and then an almost linear decrease 
of s while the neutron skin thickness of the isotopes in- 
creases. In Ref. [sl] we have studied a formation of a 
neutron skin in tin isotopes with smaller A where very 
poor experimental information is available. For instance, 
a large uncertainty is shown to exist experimentally in the 
neutron skin thickness of ^^*Sn, i.e., its value varies from 
0.1 to 0.3 fm depending on the experimental method. 
Our theoretical prediction AR — 0.17 fm for this nucleus 
is found to be within the above experimental band. A 
similar approximate linear correlation between AR and 
Po for Sn isotopes is also shown in Fig. IH^b). There- 
fore, the neutron skin thickness of these nuclei can be 
extracted once the pressure po at saturation density is 



known. The asymmetric compressibility AK given in 
Fig. in^c) is less correlated than pq with AR within the 
Sn isotopic chains. 

Figs. [7] and [5] illustrate the evolution of the symmetry 
energy with different Skyrme forces and the change of the 
neutron skin for Sn isotopes in SLy4 case. It is seen from 
Fig. [7] that the symmetry energies deduced from different 
Skyrme parametrizations vary in the interval 27-30 MeV. 
Similarly to the Ni isotopic chain (see Fig. [3]), the LNS 
force yields values of s that overestimate the correspond- 
ing values of the other three Skyrme-type interactions. 
Also, the results for the symmetry energies calculated 
with both SG2 and Sk3 forces are comparable. At the 
same time the slope of all curves remains almost constant. 
As can be seen from Fig HI the values of AR for ^^°Sn 
and ^'^^Sn obtained from our HF-I-BCS calculations with 
SLy4 force fit very well the corresponding neutron skin 
thicknesses of 0.23 ±0.04 fm and of 0.24±0.Q4 fm of these 
nuclei derived from pygmy dipole resonances (23 | . More- 
over, both experimental and theoretical values follow a 
trend established by a measurement of the stable Sn iso- 
topes ^Ag j62|]. On the other hand, the extended region 
of Sn isotopes in the present study gives a good base for 
further measurements and systematic microscopic model 
calculations in order to assess the isospin dependence of 
the nuclear EOS. 

In Fig. iniis given the evolution of the ratio AR/I for 
Sn isotopes. A clear signal that a formation of a neutron 
skin can be expected to start at A > 132 is seen in favour 
of the magic nature of ^^^Sn. Additionally, there is an in- 
dication for a weak shell closure at ^''°Sn (iV = 88). As in 
the case of Ni isotopes, the proportionality of the neutron 
skin thickness AR to the asymmetry parameter / follow- 
ing from the liquid-drop mass formula [Eq. (^T])] leads to 
a gradual linear increase of AR/I with A also within the 
Sn isotopic chain. The corresponding numerical values 
of this ratio turn out to occupy some intermediate range 
comparing to our HF-I-BCS results. 

The kinks displayed in Figs. [T][3] by the Ni isotopes 
and in Figs. |B][5] by the Sn isotopes can be attributed 
to the shell structure of these exotic nuclei. Indeed, the 
isotopic chains of Ni and Sn are of particular interest 
for nuclear structure calculations because of their pro- 
ton shell closures at Z = 28 and Z = 50, respectively. 
They also extend from the proton drip line that is found 
nearby the double-magic ^^Ni (see, for instance, the dis- 
cussion in Ref. ^2| for a possible proton skin formation) 
and ^"^Sn nuclei to the already (3 unstable neutron-rich 
double- magic ^^Ni and ^^^Sn isotopes. In the case of ^*Ni 
{N = 50 shell closure) the filling of the 139/2 orbit is com- 
pleted. Beyond this isotope the 2^5/2 subshell is being 
filled and a thick neutron skin is built up. Similar picture 
is present for the Sn isotopes. Here, one finds a sudden 
jump beyond iV = 82 where the I/111/2 shell is filled and 
the 2/7/2 subshell becomes populated. The situation in 
Pb isotopes shown in Fig. [5] is different from those of Ni 
and Sn isotopes-no kinks appear in the Pb chain consid- 
ered. In the case of ^^^Pb (N = 126) the 3pi/2 orbit is 
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FIG. 5: HF+BCS neutron skin thicknesses AR for Pb isotopes as a function of the symmetry energy s (a), pressure po (b), 
and asymmetric compressibihty AK (c) calculated with SLy4 force. 
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FIG. 6: (Color online) Same as in Fig. [T] but for Sn isotopes. 



filled and above this nucleus the first occupied level is 
2(79/2- As a result, the bulk and surface contributions to 
the neutron skin thickness of neutron-rich Sn and Pb iso- 
topes reveal an opposite effect, namely the surface part 
dominates the bulk one in tins, while for Pb isotopes the 
bulk part is larger p6l |. 

As it has been shown by Brown and Typ el (se e 
mill), and confirmed later by others [i,[2llill63,[3, 
the neutron skin thickness pO)) calculated in mean-field 
models with either nonrelativistic or relativistic effective 
interactions is very sensitive to the density dependence 
of the nuclear symmetry energy and, in particular, to 
the slope parameter L (respectively to the pressure po 
[Eq. (IH)]) at the normal nuclear saturation density. We 
note that while the basic idea behind studies of this 
type of correlations, namely between the neutron skin 
thickness and the symmetry energy (or the slope of the 
neutron EOS), is to constrain the symmetry energy in 
bulk matter from experimental results in finite nuclei, 
the main aim of our work is to extract quantitative in- 



formation about s, pq and AK values for finite nuclei in 
CDFM. Moreover, so far in our work we presented an- 
other type of correlation, namely how these quantities 
are related to AR within a given isotopic chain apart 
from the same one for a specific nucleus using different 
theoretical models. Along this line we analyze the lat- 
ter correlation and results for several Sn (^'^^Sn, ^■^^Sn, 
^^^Sn) isotopes and ^^Ni nucleus are shown in Fig. [TO] 
Using the four sets of Skyrme interaction parameters, a 
linear fit to the correlation between AR and po is per- 
formed to illustrate qualitatively the above correlation. 
It is seen from Figs. rTOT a) and rTOT d) that in the cases of 
double-magic ^"^^Sn and ''^Ni nuclei AR almost linearly 
correlates with pq. If one goes from ^^^Sn nucleus fur- 
ther to unstable nuclei within the same Sn isotopic chain 
the correlation becomes weaker, as it can be seen from 
Figs. fTOT b) and [TOT c) . being almost the same in the case 
of neighbouring ^"^^Sn nucleus and poorly expressed in 
the case of far away ^^^Sn isotope. Due to the limited 
number of Skyrme parametrizations used in our theo- 
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retical approach, the results presented in Fig. [TOl cannot 
provide a definite constraint of the slope {po value) of 
the symmetry energy. Obviously, more complete confir- 
mation of the existing correlation between the neutron 
skin thickness and the pressure for a given nucleus can 
be drawn when large set of nuclear models may be con- 
sidered. 

Finally, we would like to give estimations for the sym- 
metry energy, pressure and incompressibility of asym- 
metric nuclear matter corresponding to the Brueckner 
EOS provided by the energy-density functional reported 
in Refs. [13, El]. They are calculated at equilibrium den- 
sity Po = 0.204 fm~'^ for which the energy per nucleon Vb 
[Eq. dn])] has a minimum. The latter {E = -16.57 MeV) 
is reached from the Brueckner formula at Fermi momen- 
tum kp = 1.446 fm"i and at the value of po = 0.204 fm"^ 
used by us, as shown in Refs. |47[ Issj . This value of po is 
considered to be quite realistic [83]. The use of Eqs. 
(P7)) leads to the following estimations: S"^*^ = 35.07 
MeV, p^^^ = 1.82 MeV/fm3, and AK^^'^ = -393.85 
MeV. It is seen that the symmetry energy is few MeV 
above the experimental range 28 MeV < S"^*^ < 32 MeV 
[11] and the values deduced from our calculations for all 
finite nuclei considered. Also, the values of p^*^ and 
^j^NM gligiitly overestimate the corresponding ones ob- 
tained in finite nuclei. Actually, the problem that the 
surface symmetry energy is very poorly constrained in 
the current energy density functional parametrizations, 
as well as by available phenomenological estimates, has 
been recently indicated in Ref. ^]. To illustrate this 
fact we would like to present here new symmetry en- 
ergy values derived from the extended LDM. By taking 
Sv = 35.07 MeV which, however, comes out from the 
Brueckner EOS and a smaller value of ys = 1.1 (corre- 
sponding to the fit with LDM formula in [sol]), Eq. p2)) 
gives values for the symmetry energy of 27.89 MeV for 
^^Ni, 28.84 MeV for "2Sn and 29.58 McV for ^ospb which 
are in a better agreement with the results obtained from 
Eq. (j28p . Studying relationships between the neutron 
skin and various properties of finite nuclei and infinite nu- 
clear matter, it follows from the comparison that our mi- 
croscopic theoretical approach produces results for these 
properties consistent with those provided by the Brueck- 
ner energy-density functional which may reduce the un- 
certainties of their future experimental extraction. 



V. CONCLUSIONS 

In this work, a theoretical approach to the nuclear 
many-body problem combining the deformed HF+BCS 
method with Skyrme-type density-dependent effective in- 
teractions [tJ] and the coherent density fiuctuation model 
[49I HO] has been used to study nuclear properties of fi- 
nite nuclei. For this purpose, we examined three chains 
of neutron-rich Ni, Sn, and Pb isotopes, most of them 
being far from the stability line and representing an in- 
terest for future measurements with radioactive exotic 
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beams. Four Skyrme parametrizations were involved in 
the calculations: SG2, Sk3, SLy4 and LNS. In addition 
to the interactions used in Ref. [63, the LNS effective 
interaction was used in the present paper, because it has 
been built up from the EOS of nuclear matter including 
effects of three-body forces [78]. 

For the first time, we have demonstrated the capabil- 
ity of CDFM to be applied as an alternative way to make 
a transition from the properties of nuclear matter to the 
properties of finite nuclei investigating the nuclear sym- 
metry energy s, the neutron pressure po and the asym- 
metric compressibility AK in finite nuclei. This has been 
carried out on the base of the Brueckner energy- density 
functional for infinite nuclear matter. Mainly, the sym- 
metry energy has been determined from analyzing nu- 
clear masses within liquid-drop models, while much more 
effort has been recently devoted to extracting the value of 
its slope parameter. Instead of this, in the present work 
we applied the CDFM scheme. One of the advantages 
of the CDFM is the possibility to obtain transparent re- 
lations for the intrinsic EOS quantities analytically by 
means of a convenient approach to the weight function. 
The key element of this model is the choice of density 



distributions which were taken from self-consistent de- 
formed HF-I-BCS calculations. 

We would like to emphasize that our new method al- 
lows one to estimate the symmetry energy and, there- 
fore, the symmetry-energy coefficient in finite nuclei thus 
avoiding the problems related to fitting the Hartree-Fock 
energies to a LDM parametrization. 

We have found that there exists an approximate linear 
correlation between the neutron skin thickness of even- 
even nuclei from the Ni {A = 74-84), Sn {A = 124-152), 
and Pb {A = 206 — 214) isotopic chains and their nuclear 
symmetry energies. Our HF-I-BCS calculations lead to 
a symmetry energy in the range of 27-30 MeV, which is 
in agreement with the empirical value of 30 ± 4 MeV 
[S^ l and with the results of Nikolov et al. [s^. In 
the cases of Ni and Sn isotopes the symmetry energy 
is found to increase almost linearly with the mass num- 
ber A till the double- magic nuclei (^®Ni and ^^^Sn) and 
then to decrease, being larger when LNS force is used. 
A similar linear correlation between AR and po is also 
found to exist, while the relation between AR and AK 
is less pronounced. The same behavior containing an 
inflexion point transition at specific shell closure is ob- 
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served for these correlations. The calculated values of 
Po = 1.36 — 1.68 MeV/fm"^ obtamed in our calculations 
lead to values of the slope parameter L = 26 — 32 MeV 
that are in agreement with other theoretical predictions 
(see, for example @|). Hence, the study of other two 
EOS parameters (whose determination is more uncertain 
than the symmetry energy), namely the slope and the 
curvature of s, performed in our work may provide more 
constraints on both the density dependence of the nu- 
clear symmetry energy and the thickness of the neutron 
skin in neutron-rich medium and heavy nuclei. 

In the present work, the role of the relative neutron- 
proton asymmetry / on the size of the neutron skin AR 
was studied for Ni and Sn isotopes. For this purpose the 
results for the ratio AR/I were presented and a compar- 
ison with its estimation using the extended liquid-drop 
mass formula [Eq. (|4T|) ] was made. The HF-I-BCS re- 
sults show that a formation of a neutron skin can be 
expected to start at ^ > 78 and A > 132 for Ni and Sn 
isotopes, respectively. Although the proportionality that 
exists between AR and / in (|^T|) produces smooth change 
of AR with A, the size of the neutron skin in both cases 
is consistent with the results of other calculations with 
the SLy4 parametrization [63|, as well as with results 
of spherical self-consistent HFB calculations with finite- 
range forces of the Gogny type . It is confirmed that 
the neutron skins in nuclei with a large neutron excess 
near and beyond the drip line are very clearly related to 
the asymmetry parameter /. However, the shell effects, 
which are always present in mean-field calculations, pro- 



duce deviations from the linear dependence between AR 
and / in these nuclei. 

Summarizing, in this work we use a microscopic theo- 
retical approach to study important macroscopic nuclear 
matter quantities in finite nuclei and their relation to 
surface properties of neutron-rich exotic nuclei. We have 
shown that this approach gives a very reasonable descrip- 
tion of the neutron skin thicknesses of several Sn isotopes 
which have been experimentally measured. A good agree- 
ment is achieved with other theoretical predictions for 
the volume symmetry energies of Ni isotopes, as well as 
for the pressure and asymmetric compressibility of nu- 
clei within the Pb isotopic chain. The capability of the 
present method can be further demonstrated by taking 
into consideration Skyrme-type and relativistic nuclear 
energy-density functionals. In this way, more systematic 
analysis of the relationship between the symmetry en- 
ergy and the neutron skin thickness will be carried out 
and more definite conclusions on its isotopic sensitivity 
will be drawn. 
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